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Abstract 



We prove the existence of asymptotic two-soliton states in the Fermi-Pasta-Ulam 
model with general interaction potential. That is, we exhibit solutions whose difference 
^ I in £^ from the linear superposition of two solitary waves goes to zero as time goes to 

infinity. 

1 Introduction 



In a series of recent works, (e.g. [TJ [2l [3l [4]) Martel, Merle, and collaborators have 
studied solutions of the generalized Korteweg DeVries (gKdV) equation 



ut = iu,, + uPU xeR; p = 2,3,4,5. (1) 

I It is a long standing observation that gKdV admits a family of solitary waves u{x, t) = 

QP ■ Uc{x — ct) in which larger waves move faster than smaller waves. When p = 2, 3 the 

system is completely integrable and one may construct 



. ^ _ pure n-soliton solutions^ i.e. solutions whose difference from the linear su- 

■ perposition of n independent solitary waves, each moving with its own char- 

H \ acteristic speed, goes to zero in e.g. as time goes to either forward or 

backward infinity. 

An important and long-standing observation about these solutions is that when a larger, 
faster soliton overtakes a smaller, slower soliton the collision is clastic - there is no 
energy lost to dispersion. It has long been an open cj^uestion as to what extent this 
phenomenon persists for non-integrable Hamiltonian systems. When p = 4, 5 the gKdV 
equation is believed to be non-integrable and thus the inverse scattering transform and 
other methods for constructing pure n-soliton solutions are unavailable. However, one 
may (and Martel does [3]) construct 

asymptotic n-soliton states, i.e. solutions whose long time behavior (in for- 
ward but not necessarily backward time) is precisely the same as that of the 
linear superposition of n independent solitary waves, each moving with its 
own characteristic speed. 

The construction of such solutions does not address the question of whether or not the 
collision is elastic, but it is a crucial step in the construction of 

nearly n-soliton solutions: Solutions whose long time behavior is precisely 
the same the the linear superposition of one collection of n solitary waves 
as time goes backward to infinity and precisely the same (when measured in 
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a norm which ignores dispersion) as the hnear superposition of a different, 
and hkely lower energy, coUection of sohtary waves as time goes forward to 
infinity. 

The difference in the energy of the two asymptotic 71-sohton states accounts for the 
energy which has been lost to dispersion as the solitons which comprise the nearly 
n-soliton solutions collide inelastically. 

One limitation of the gKdV equation is that it supports wave propagation in only 
one direction. More realistic models, like the equations for water waves allow waves to 
propagate in any direction. Another class of equations which supports two-wave wave 
propagation is the Fermi-Pasta-Ulam system of Hamiltonian differential equations: 

f (S* - I)p 

(2) 

p={I-S~^)V'{r) 

where u = {r,p) £ x fi, {Sx)n — Xn+i is the shift on i"^, and F : R — > R is an 
interaction potential which satisfies V{0) = V'{0) = 0, ^"(0) > 0, and V"'(0) ^ 0. (In 
fact, it is simple to show that one can assume, without loss of generality that V"{0) = 1 
and we will assume this normalization from now on.) An extensive theory of small 
amplitude, long wavelength solitary waves has been developed for ([2]). In particular, 
waves with amplitude ^ e^, wavelength ~ and wavespeed ~ 1 + exist [5], are 
exponentially localized and supersonic [5] , and stable in a variety of senses of the word 
[21 [8]. We shall denote the solitary wave with speed c = ±(1 + e^/24) by Uc = {rc,Pc)- 

Remark. Of course, any translate of the solitary wave Uc is also a solitary wave of the 
same speed. For definiteness, we assume that Uc is the member of this family for which 
the profile rc has its (unique) maximum at zero. 

We have recently investigated the long-time behavior of counterpropagating solitary 
waves in the small-amplitude long- wavelength regime of the FPU model [9j . In particu- 
lar we have demonstrated that solutions exist which are close to the linear superposition 
of two solitary waves for all time and that moreover these solutions are asymptotically 
stable with respect to perturbations which arc exponentially localized and orbitally sta- 
ble with repsect to perturbations which live in P . The purpose of this note is to prove 
the existence of asymptotic counterpropagating 2-soliton states in the long-wavelength, 
low-amplitude regime of the FPU model. 

Theorem 1.1. The FPU model admits asymptotic two soliton states. That is, given 
c_ sufficiently close to and less than —1 and c+ sufficiently close to and greater than 1, 
there is an initial datum uq G £^ such that the solution of ^ with initial data u(0) = uq 
satisfies 

lim ||u(t)-Ue+(--c+i-7+(i))-Uc_(--c_t-7_(i))|| =0 (3) 

t^oo 

for some functions 7± which satisfy limt^oo 7(t) = 0. 

Remark. Here, and in what follows, \\ ■ \\ refers to the £^(Z) norm. 

At the heart of this existence result is the following orbital stability result which 
was proven in [9j. 

Theorem 1.2. There are positive constants C, Tq, b, 6q, and Eq such that for any 
e G (0,eo) the following hold. Suppose that 

\\uo{-) - Wc+(- - T+) - Uc_{- - t1)\\ <S <6o 

with t1 < —T < T < and T > Tq. Then there are smooth real valued functions of a 
real variable C-|_, c_, r-|_, and t_ such that the solution u of ^ with u(0) — uq satisfies 

\\u{t, ■) - «,^(,)(. - r+(t)) - u,_^t}{- ~ r.m < + Ce""^^ (4) 
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(5) 
(6) 



for all t > Moreover, 

\c+{t) - c+{t„)\ + |c_(t) - c_(to)| < Ce-^S + Ce-"'^. 
Finally, i/ie /imiHimt_,oo(c+(t), c_(t), f-|_(t), f_(t)) exists with 
\T±{t) - c±it)\ < Ce-^/^S + Ce~^^^ 

holding uniformly for t > and limt_»oo |7'±(0 ~ c±(i)| =0 

Theorem 11.21 follows by combining the results of Theorem 2.2, Theorem 4.5, and 
Lemma 3.7 from [9]- 

2 Preliminaries 

In this section we collect several lemmata which streamline the proof of Theorem 11.11 
The first such lemma provides an estimate on the residual of the linear superposition of 
two solitary waves which is valid on exponentially long time scales. The application of 
this lemma is a key step in the proof of Theorem ll.il The other lemmas in this section 
serve to verify the smallness condition ([8]) in the hypotheses of Lemma 12.11 



Lemma 2.1. Consider the discrete wave equation which is driven, from rest, by small 
inhomogeneous and nonlinear terms: 



' p=iS-I)TT 

7r= {I - S-'^)p + g{TT,p,t) 
[ p{0) = 7r(0) = 



(7) 



Let N and c± be given real numbers. Define N± := c±N. Let {x,y) :~ '^jtzz^jUj 
denote the usual P' inner product and define weighted inner products by {x,y)± := 



ergy" of the solution (p, tt) of ^ by 



"Cjyj. Define the "en- 



m 2(l|pWf + llpWf+ + llpWII- + lk(t)f 



Mt)\\\ + Mm-) 



Given positive constants C , C , rj, and b, there is a positive constant Eq such that as 
long as g satisfies the smallness condition 



J2 \\gi7:,P,t)UHt)\\o < C [(£3/2 + E(i)^/')E(0 +e^+''e-''^(^-*)' 

«G{^2, + ,-} 

for all t < N and for some e S (0,eo), then the energy E satisfies the estimate 

E(t) < C'e^e-''^^^-*' ie[0,iV] 

Proof. Compute 
E{t) 



(8) 



(9) 



C+£3/2( 


pmi + 


||7r(t)||^)-c_£3/2(||p(t)l|2_4 






- /)7r) + 


{n,{L-S-')p) 


(a 


+{p, (s 


-/)7r)+H 


-{n,{I-S-')p) + 


[Hi 


MpAs 


-/)7r)_H 




{iv 






+ (^,.9)- 


{v 



Note that (i) is bounded by (c+ + c_)£'^/^E(<) and (m) is identically zero. Since the 
shift S is not self-adjoint on the weighted spaces, we estimate more carefully 

|(m)| = - S-')p,7:)+\ = \ie-'^'''- 1){S-' p,7:)+\ < Coe'/^pU\\nU 

and similarly < Coe^/^\\p\\^\\7r\\^. The estimate |(i;)| < Ce^/Se + Ce'*+''e-''^l^^*l 
follows immediately from ([8]) and the Cauchy- Schwartz inequality as long as E{t) < Ce^. 
Choose and fix (7 > and define T := sup{i e [0, N] \ E(i) < Ce'^] so that 

E < Cse^/^E + C'e4+,,g-6e(Ar-t) ^ ^ JQ^ 

E(0) = 

Apply Gronwall's inequality to see that 

< _£0^e^^-'''*-^^^(e(''^-^-^''')* - 1) iG[0,T] 

where in the last inequality we have taken e sufficiently small. Thus a posteriori we 
may take T ^ N. This establishes © and thus completes the proof. □ 

The following simple lemma establishes that the weighted space with norm E embeds 
compactly into 

Lemma 2.2. Let a > be given. Let X CZ i'^ en-joy the estimate 
R{X) sup x]{l + e"^ + e""^) < oo. 

Then X is compact in 

Proof. We show that for each 77 > there is an iV = N{ri) and a discrete set D ^ D{ri) = 
{xi, • ■ • xjv} such that the 77-neighborhood of D contains X. Let 77 > be given. Choose 
M so large so that for any x G X, we have X]|j|>M I-^j I^ ^ '7- Such a choice is possible 
because R{X) is finite. Now let D be a discrete set such that the jy-neighborhood about 
D contains the ball of radius R{X) in {x \ = for \j\ > M} ^ M^flz+i g^^^,;^ 
a set exists because R^*^+^ is locally compact. Note that the 77-neighborhood of D 
contains X as desired. □ 

The following lemma establishes that a key cross term is exponentially small. We 
make use of this lemma in verifying the smallness condition ([8]) in the hypotheses of 
Lemma 12. f I 

Lemma 2.3. Let yj{t) r^^iJ ~ c+{N ~ t)y^_{j - c^{N - 1) - a) for some a £ (0, f). 
Then there is ah > Q which may he chosen independent of e and a such that the estimate 

\\y[t)U + \\y[t)\\^ + \\v{t)\\<Ce''^e-'<''-'^ t<N (10) 

is valid. 

Proof. It is known from j6j that the wave profile and its derivative enjoy the 
exponential estimates 



-efco|?| 
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for some positive constants C and feg which may be chosen independently of e. Thus 

\yj{t)\ < Ce5g--f,oeb-c+(JV-t)|g-bo6|j-c-{A'-t)| 

There are three regimes to consider: (i) j > c+{N — t), (ii) c_ {N — t) < j < c+{N — t) 
and (iii) j < c_ {N — t) . We first compute the i'^ norm: 

]<c-{N-t) c-{N-t)<]<c+{N-t) j>c+{N-t) 



fe>0 



E 

C-(W-t)<J<c+(JV-t) 



-2boel 



l>0 



After the first inequality we have changed indices k = {N — t)—j and I = j — {c-^{N — 
t)). In the last hne we have used the fact that sup^>g xe^'^^ < Ca^^ for any a > and 



in particular for a = boe. We have also used the fact that for a > 0, X]fc>o ^ 
with a = 2bQe. Wc now compute the || • ||+ norm: 



gl0g-2e3/2c+(jV-t)g-2hoe[(c+-c_)(W-t)] 



fe>0 ' 



„2e3/2 (fe+c+ (JV-t)) p-4fcoefc 



^^^^ g2£^/='(c_(7V-t)-0g-4ho£i _|_ ^ 



je(c_(iV-t),c+(Ar-t)) ' 



< C'eiOe-^[(^+-=-)(^-*)l + [(c+ - c_)(7V - t)]) 

Here we have used the fact that e'^/^c+ < ^(c+ — c_) for small e to ignore all of the 
powers of e'^/^ and essentially reduce this computation to that of the £^ norm above. 
The estimate for is similar. After taking square roots, we obtain (fTO|) : this 

concludes the proof. 

□ 

The following lemma will be used to establish that the inhomogcneous terms which 
drive the evolution of the perturbation of two well-separated solitary waves are bounded 
by a cross term of the form shown to exponentially small in Lemma 12.31 We leave the 
proof as an exercise. 

Lemma 2.4. Let f : R ^ R be and fix 0. Then the function G : _ 



G(a;,y) := < 



r f{x+y)-f{x}-f{v) 

xy 

/'fa) 
V 



X 

I /"(o) 



x,y 

x = Q,y^O 

x^Q,y = 
X = y = 



given by 



(11) 



is Lipschitz continuous in any bounded set in R^. 



3 The proof of Theorem 11.11 

We are now ready to prove Theorem ll.il 

Proof of Theorem \l.l\ In broad strokes, the proof proceeds as follows. Ultimately we 
wish to have a pair of solitary waves u^^ and Uc_ moving away from each other to the 
right and left respectively. However, for the first step in the proof we consider their 
evolution toward each other from a widely separated initial state. Note that because of 
the time reversibility of the FPU equations, if the solution Wc+ (■, t) — {rc^ (•, i), Pc+ (•, t)) 
is a solitary wave which propagates to the right with speed c+ > 1, then U^^{-,t) = 
(rc^ (•, — t), — Pc+ (■, —t)) is a solitary wave propagating to the left with speed — c+. Like- 
wise we introduce a solution which is the "time reversal" of the solitary wave Uc_ 
and represents a solitary wave propagating to the right. With this notation we consider 
a sequence of solutions of the FPU equations, {u^{-,t)}, with initial conditions 

M^(-,0) = U"{- - c+N,0) + U^_{- - c^N,0) . 

Let this initial data evolve for time N — Tq, at which point it may be written as the linear 
superposition of colliding waves which remain separated by a distance of (c+ — c_)To 
plus an error term. We prove that this error term is uniformly small in A'^ in a space 
which embeds compactly into f^. Thus the sequence u^{N — Tq) converges in i'^ along a 
subsequence Nj to some limit. Use reversibility to make the waves in this limiting initial 
datum move apart rather than together and call the result uq; its forward evolution u 
is the solution that we seek. At this point we rely on Theorem 11.21 which says that so 
long as initial data starts close to the sum of two well-separated solitary waves which 
are moving apart from each other, the evolution remains close to the sum of these two 
solitary waves for all time. Moreover, the distance to the sum of two solitary waves is a 
constant multiple of the initial distance plus a term which is exponentially small in the 
separation. 

We now use the fact that the solutions u^^ {t) converge to the solution u{t) whose 
initial condition is uq. Since u^^ (t) is equal to a pair of solitary waves when t = Nj, we 
find a sequence of times tj, tending toward infinity, such that for t > tj, u{t) is more 
and more closely approximated by a pair of solitary waves. Allowing t to go to infinity 
we find that u{t) converges to a pair of solitary waves with no radiation. We note that 
this general method of using apriori bounds to derive solutions with specified properties 
at t — > oo has been used in a dissipative context in Hamel and Nadirashvili |10| , and in 
the gKdV equation by Martel [5] . 

We now present the details. Let u^{j, 0) := U^^{j - c+N) + {j - c_N) and let 
u^{-,t) be the evolution of u^{-,0) under the FPU flow. Define the residual (p^,7r^) 
by u^{j, t) = C/« (j - c+N + c+t) + Ut U - c-N + c^t) + t), ^^(j, t)). 

Suppress the superscript N . The residual (p, tt) satisfies the equation 

' p={S-I)^ 

< ,V = (/ - 5-1) [V'ir,^ + r,_ + p) - F'(r,J + V'{r,_ )] 
^ p(0) = ^(0) = 

where Uc = (tcPc) and rc_^_ and rc_ are evaluated at j — c+N + c+t and j — c^N + c^t 
respectively . 

Rewrite the right hand side of the equation for tt as 

(1 - S-')p + (1 - S-'Wir,^ + r,_ +p)- r(r,J - V'ir,_)~ p). 



fi 



The quantity V'{rc_^_ + rc_ + p) — V'{rc^) — V'{rc_ ) — p will play the role of g{n, p, t) in 
Lemma l2.ll and we begin by decomposing it as: 

V'ir,^+r,_+p)-V'{r,^)^V'ir,_)-p = V {r,^ + r,_ + p) - V {r,^ + r,_) - V" (r,^ + r,_)p (i) 

+ iV"{rc^+rc_)-l)p iii) 
+V' (r,^ + r,„ ) - y (r,^ ) - V (r,_ ) {in) 

From Taylor's theorem with remainder we sec that < Cp^. So 

\\ii-s-'m\u<c\\px (12) 

for o G In like fashion, we observe that V"{0) = 1 that |(m)| < C\{rc^+rc_)p\ 

and hence 

11(1 < +r,Jp|U < C(||r,J|,^ + \\r,_\U^)\\p\U < Ce^p\U (13) 

The estimates on term (iii) are somewhat more subtle. First use Lemma [2.4l to conclude 

{I - S-^){tti) < C{I ~ S-^)r,^n_. 

Next, note that 

(/ - 5-i)r,^ (■ - e+)r,_ (. - C_) = r,^ (• - e+)(/ - S-')r,_ )i- - C-) + (^^ (• - - S'^M - 

= (• - e+)rL {■-(-) + (• - - 1)^:;+ (• - c+). 

Here ^± = c±(A^ — t) and C± G (C± ^ 1; ^±); in the second line we have used the Mean 
Value Theorem. It now follows from Lemma [2731 that 



II (/ - 5-i)r,,r,_ lU < Ce^/^e-''^^^-*), t < N, 

for any o G {£^, +, -}. 

Multiplying the above bounds by ||7r||o and then summing over o £ +, — }, we see 
that dS]) is satisfied and thus if we fix some Tq > 0, (p^(A^-To), 7r^(A^-ro)) is bounded 
uniformly in the space i"^ n i\ n ("i. We established in Lemma that i'^ n £\ n it 
embeds compactly in it follows that (p^(iV - To), 7r^(iV - Tq)) converges strongly 
in along a subsequence Nj to some limit (p, tt). 

Note that the limiting solution consists of a pair of solitary waves at positions c+Tq 
and c_To, plus a small remainder. We will assume that Tq is large enough that the 
"overlap" of the tails of the waves remains small. As described above we now consider 
this limiting solution as the initial condition for a solution in which the two solitary 
waves move apart from one another. Once again, we make use of the reversibility of the 
FPU equation (and the definition of C/^?^ and U^_) and write 

uoU) = Uc^{j - c+To) + Uc_{j - c_ro) + {p{i), -7r(j)) 

and let u be the solution of ^ with u(-,0) = uq. Note that due to reversibility and 
continuous dependence on initial conditions u{-,t) is close to u^{-,N ~Tq ~ t), where 
u^^ represents the "time reversed" analogue of u^^ - i.e. we insert a minus sign in front 
of second component of the solution (the p-component) to account for the reversal of 
time. In particular, at time t, u will be close to a pair of "outward" propagating solitary 
waves located at c± {Tq + t) . 

The argument proceeds by using the orbital stability result. Theorem 1 1.2 1 to control 
the forward evolution of u for all time. Due to the presence of the e^*"^^ term on the 
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right hand sides of ([4]) and ([5]), it is not useful to regard uq as the initial condition for 
u. Instead, we consider a sequence of times Tfc — > oo and regard u{-,Tk) as the initial 
condition for u. The exponentially small terms in ^ and ([5]) arc now order and 
upon letting k ^ oo they vanish. The remainder of the proof makes this argument 
precise. 

Choose sequences {Sk} and {T^} such that 5^ — > and Tk ^ oo. Then by continuous 
dependence on initial conditions there exists N^. > T^. such that 

\\ui;Tk)-u'''i;Nk-Tk~To)\\ <6k (14) 

Now write 

u^'-i; Nk ~Tk)^ uc+ (• + c+ {Nk -Tk- To)) + Uc_{- + C- {Nk - Tk - To)) 
+ (•, Nk - Tk - To), -7^""^ {■,Nk -Tk- To)) 

By Lemma [2. II we have 

IKp^n--^^ - Tk - To), -7r^H-,A^fe -Tk- To))\\ < Ce^e-''^^''+^"^ (15) 

Thus 

\\ui-,Tk)-Uc^i-+c+{Nk-Tk)-c+Nk)+Uc_{-+c-{Nk-Tk)-c-Nk)\\ < Sk+Ce^'e-'"^'' 4 

(16) 

Now apply Theorem 11.21 with uo = u{-,Tk). Then for all t > Tk we have functions 
Cj. and t![ such that 

Mt) - t*,. (,)(. - T\{t)) - u,.j^^{- T^Lm\ < Ce-^/Hk + Ce-^^^ (17) 

holds. It remains to replace the /c-dependent modulation parameters r^i and with 
fc-indepcndent quantities. 

In light of ([5]) and the fact that c i-^ is smooth, we may replace c\{t) with 
c± at the cost of an additional Ce~'^5k + Ce~^'^'' on the right hand side. Now define 
7i(t) := r^(i) - c±t. In light of © and ® we see that |7|(t)| < Ce'^Sk + Ce"^^'^''. 
In particular, 7^ ^ uniformly as /c — > 00. Taking A: — > 00 we find that u converges to 
a pair of solitary waves. This completes the proof. □ 
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